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Dust formation in the winds of AGBs: the contribution at 
low metallicities 
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ABSTRACT 

We present new models for the evolution of stars with mass in the range IMq ^ M ^ 
7.5M0, followed from the pre-main-sequence through the asymptotic giant branch 
phase, until most of their envelope is lost via stellar winds. The metallicity adopted is 
Z = ?,x 10~^ (which, with an a— enhancement of +0.4, corresponds to [Fe/H] = —2). 
Dust formation is described by following the growth of dust grains of various types as 
the wind expands from the stellar surface. 

Models with mass M ^ ZMq experience Hot Bottom Burning, thus maintaining 
the surface C/0 below unity. Unlike higher Z models, the scarcity of silicon available 
in the envelope prevents the formation of silicates in meaningful quantities, sufficient 
to trigger the acceleration of the wind via radiation pressure on the dust grains formed. 
No silicate formation occurs below a threshold metallicity oi Z — 10^"^. 

Low-mass stars, with M ^ 2.5M0 become carbon stars, forming solid carbon dust 
in their surroundings. The total dust mass formed depends on the uncertain extent 
of the inwards penetration of the convective envelope during the Third Dredge-Up 
episodes following the Thermal Pulses. However, provided that a minimum abundance 
of carbon oiX{C) '^ 5 x 10"'^ is reached in the envelope, the results turn out to be fairly 
independent of the parameters used. Carbon grains have sizes 0.08 mum < ac < 0.12 
/im and the total amount of dust formed (increasing with the mass of the star) is 
Mc^ (2-6) X 10-4 Mq. 

Our results imply that AGB stars with Z = 3x 10"^ can only contribute to carbon 
dust enrichment of the interstellar medium on relatively long timescales, > 300 Myr, 
comparable to the evolutionary time of a 3Mq star. At lower metallicities the scarcity 
of silicon available and the presence of Hot Bottom Burning even in M< 2 Mq, prevents 
the formation of silicate and carbon grains. We extrapolate our conclusion to more 
metal-poor environments, and deduce that at Z < 10~* dust enrichment is mostly 
due to metal condensation in supernova ejecta. 

Key words: Stars: abundances - Stars: AGB and post- AGB. ISM: abundances, dust 



1 INTRODUCTION 

The importance of intermediate mass stars on the asymp- 
totic giant branch (AGB) and super asymptotic giant 
branch (super-AGB) as source of cosmic dust is cur- 
rently a hig hly debated topic. R e cent c hemical evolution- 
ary models (jValiante et al.l 12003 . |201l|) have shown that 
AGB stars give important contributions to dust enrich- 
ment even at early cosmic epochs, when traditionally core- 



collapse supernovae (CCSNe) have be en proposed as the 
main source of dust ([Todini fc Ferrara 2001; Nozawa et al.l 
l2003l ; iBianchi fc SchneideJ |2007| ; iDweklTioOTJ) . The recent 
detection of about (0.4 — 0.7)Mq of cold dust in the su- 
perno va 1987A by the Herschel satellite (jMatsuura et al.l 
[2OIII) confirms the results expected by theoretical models 
of dust condensation in supernova ejecta. Yet, only a rel- 
atively small fraction of the newly formed dust, probably 
less than 20%, is able to survive the passage of the re- 
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verse shock and be injected in the surrounding interstellar 
medium (ISM), e ffectively contributing to dust enrichment 
of th e host galaxy (JBianchi &: Schneideill2007l : lNozawa et alJ 
|2007|). Indeed, Spitzer observations of young supernova rem- 
nants have provided unabiguous evidence of dust formation 
and survival in the ejecta, with typi cal masses in the range 
(1 - 5) X 1O"^M0 (|Rho et al.ll2009l ). These figures can be 
considered as an indication of the effective dust yields of 
massive stars. 

At lower stellar masses (< 8Mq) grain condensation 
occurs during the AGB and super-AGB phases and the re- 
sulting dust mass depends on the progenitor stellar mass 
and metallicity through complex processes which affect the 
chemical composition and physical conditions in the stellar 
atmospheres. 

Previous investigations, based on synthetic stellar mod- 
els, have provided grids of AGB dust yields predicting dust 
masses in the range (10~* — 10~^)Mq (Ferrarotti & Gail 
2006). With the aim to make further progresses follow- 
ing the proper stellar evolutionary sequences, we have re- 
cently computed the dust masses released by stars with 
masses in the range IMq ^ M ^ SM© and metallicities 
of Z = 0.001 (5 X IQ-^Zq) and Z = 0.008 (0.4^©). We find 
that the both the total mass of dust released and its chemical 
composition are a strong function of the initial stellar mass 
and metallicity (Ventura et al. 2012a, 2012b) and that, in 
the most optimistic cases, the typical dust masses released 
by AGB stars can reach up to 7 x 1O~^M0. While the rela- 
tive importance of AGB stars and SN as sources of dust will 
ultimately depend on the star formation history and stellar 
initial mass function (jV aliante et al.ll2009l . l201ll ). in order 
to understand the contribution of AGB stars in early cos- 
mic dust enrichment we need to extend our previous models 
to lower initial metallicities. 

In this paper, we calculate a new sequence of models 
with metallicity Z = 0.0003 (1.5 x lO'^Z©), evolved through 
the whole AGB phase. Our results imply that at these low 
initial metallicities, the production of dust is strongly sup- 
pressed, suggesting that in cosmic environments with metal- 
licities Z < 5 X 10~^Zq dust enrichment is entirely domi- 
nated by supernovae. 



2 THE MODEL 

The thermal pulses (TP) phase is experienced by stars of in- 
termediate mass (M < {8 — 9)Mp)), a fter helium is consumed 
in the core (see e.g. iHerwia (120051 ) for an exhaustive de- 
scription of these phases). Helium and hydrogen burning oc- 
cur in two shells above the degenerate core. This evolution- 
ary phase is commonly known as Asymptotic Giant Branch 
(AGB). In the most massive models, whose core mass ex- 
ceeds ~ I.IA/q, carbon ignition occurs in an off-center, de- 
generate region. This carbon pulse favors the development 
of a convective flame, that moves inwards to reach the stellar 
center, leadin g to the formation of a core composed of oxy- 
gen a nd neon l|Garcia-Berro et al.ll 19971 : iRitossa et al.lll99a . 
Il999r ). These stars also experience a series of thermal pulses, 
duri ng the so ca l led Su per Asymptotic Giant Branch evolu- 
tion l|Siesdl2OOl[2O10l ). 

The surface chemistry of these stars may be altered by 
two mechanisms: Hot Bottom Burning (HBB) and Third 



Dredge-Up (TDU). HBB is active in massive AGBs, and 
consists in the ignition of proton-capture nucleosynthesis at 
the bottom of the convective envelope; it is activated when 
the temperature of this region reaches ~ 5 x 10^ K. TDU 
consists in the inwards penetration of the envelope after each 
TP; the bottom of the surface convection zone may reach 
regions processed by 3a nucleosynthesis, enriched in carbon 
and oxygen. 

The modification of the surface chemistry produced by 
HBB and TDU are different. TDU is accompanied by the in- 
crease in the surface abundance of carbon, and, to a smaller 
extent, of oxygen. The effects of HBB refiect proton-capture 
nucleosynthesis, and thus depend on the temperature at 
which HBB occurs. For temperatures below ~ 8 x lO'^K the 
main effect is the production of nitrogen at the expenses of 
carbon, whereas at higher temperatures oxygen destruction 
also occurs. 

Due to the large mass loss rates experienced, and the 
cool temperatures in their expanded, low-density envelopes, 
the winds of AGBs are a favorable site for dust production, 
via condensati on of gas molecules into s oUd grains. The in- 
vestigations bv lFerrarotti fc Gaij (|200ll . 12002 . 120061 ) set the 
theoretical framework to model dust formation in expand- 
ing AGB winds, and showed how the formation of carbona- 
ceous species (essentially solid carbon and SiC) and silicates 
(olivine, pyroxene, quartz) depends on the surface chemistry 
of AGBs, particularly on the C/O ratio. These works are 
based on a synthetic description of the AGB phase, where 
core mass, temperature at the bottom of the convective zone, 
extension of TDU, are all imposed a priori, based on the ob- 
servational calibration. This approach has the clear advan- 
tage of allowing one to follow the whole AGB evolution with 
small computational efforts; the intrinsic limit of this de- 
scription is that HBB cannot be described self-consistently. 

In our previous investigations (Ventura et al. 2012a,b, 
hereafter paper I and paper II, respectively) we presented 
new models where the AGB evolution was calculated by 
means of the integration of the whole stellar structure. We 
also extended the calculations to the super-AGB regime, ig- 
nored in previous works. Our analysis showed that the kind 
of dust formed around AGBs depends on the initial mass of 
the star: models with M J5 3Mq, whose surface chemistry 
is dominated by HBB, produce silicates, whereas their lower 
mass counterparts produce carbon-type dust. These investi- 
gations also showed an interesting trend with metallicity: in 
higher-^ models the production of silicates is favored, not 
only by the larger amount of silicon available in the enve- 
lope, but also because the HBB experienced is weaker, thus 
inhibiting the possibility of surface oxygen destruction. 

In what follows we further extend the model to lower 
initial metallicities and discuss the implications for the mass 
and composition of dust yields by intermediate mass stars. 

Dust grains are assumed to form via condensation of gas 
molecules in the wind of AGBs. The description of this pro- 
cess requires knowledge of the evolution of the main physical 
and chemical properties of the central object, and a model 
for the thermodynamical and chemical structure of the wind. 
We provide here a brief description of the st ellar and wind 
modeling, addressing the interested reader to I Ventura et al.l 
((2012a|) for a thorough presentation of the model. 
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2.1 Stellar evolution modelling 

The stellar models presented in this work have an initial 
metallicity Z = 3 x 10~^. The mixture is assumed to be 
a— enhanced, with [a/_Fe] = +0.4. This corresponds to an 
iron content [Fe/ H] — —2. The evolutionary sequences were 
calculated by means of the ATON stellar evo lution code, in 
the version described in IVentura et al.l (|l998l ). 

The range of masses involved is limited to M ^ 7.5Mq, 
because more massive objects undergo core-collapse SN type 
II explosion, thus not experiencing any AGB phasqj. The 
evolutionary models were followed from the early pre main 
sequence phase, until the phase where almost all the external 
envelope is consumed. Models with masses below 2.5Mq de- 
velop a degenerate helium core while ascending the red giant 
branch, and undergo the helium flash. The evolution of these 
models was thus re-started from the horizontal branch, from 
a post-flash stage, before the beginning of the quiescent core 
helium burning. 

The convective instability is modeled according to the 
Fu ll Spectrum of Turbulence (hereinafter FST) treatment 
bv lCanuto fc Mazzitellil l|l99ll ). In regions unstable to con- 
vection, mixing of chemicals and nuclear burning are cou- 
pled via a diffusion -like equation , following the scheme pre- 
sented in lCloutman fc Eolj (|l976l ). The overshoot of convec- 
tive eddies within radiatively stable regions is described via 
an exponential decay of velocities from the formal border of 
convection; this is determined via the classic Schwarzschild 
criterium, stating that the zone of neutral stability, where 
buoyancy vanishes, corresponds to the condition Vrad = Vad 
(where Vrad and Vad are, respectively the logarithmic tem- 
perature gradients when the energy is entirely transported 
by radiation, or for an adiabatic transformation). The e- 
folding distance I of the exponential decay of convective ve- 
locities is assumed to be i = (" x Hp, where Hp is the pressure 
scale height, Hp = P/{gp) and (^ is a measure of the exten- 
si on of the overshoot re gion. Following the calibration given 
in I Ventura et al.l (11998), we adopt (^ — 0.02 in all the evolu- 
tionary phases preceding the thermal pulses phase. The pos- 
sible extension of the extra mixing during the AGB phase is 
still an op en issue, largely debated in the literature (jKarakasI 
I2OIII : [Sta ncliffc & Lattanzio 2011). To limit the number of 
free parameters, here we only allow some extra mixing from 
the borders of the convective shell that forms during the 
thermal pulse, and discuss how the results obtained depend 
on the assumed value of (". 

Ma s s loss was modeled following the treatment by 
iBlockerl l| 19951 ). 

The majority of the nuclear reaction rates were taken 
from the NACRE compilation JAngulo et all Il999l'). The 
only exceptions are the 3a cross sect ions (J Fvnbo et al. 
[20^, the rate of the ^^C + a reaction l|Kunz et al.ii200al . 
and some proton-captures, not relevant for the present in- 
vestigation. 



^ This limit is indeed partly dependent on the assumption of 
some extra mixing from the border of the convective core during 
the Main Sequence phase. If overshooting is neglected, the highest 
mass not undergoing SNII explosion is ~ 9.5Mq 



2.2 Dust production and the stellar wind 

The scheme adopted to model the structure of the wind at 
a given phase during the AGB evolution is the same as in 
paper I, ba sed on the de scription given in the series of pa- 
pers bv iFerrarotti fc Gail (2001, 2002, 2006). Here we briefly 
recall the schematization adopted. 

The wind is assumed to expand isotropically from the 
stellar surface. The physical quantities characterizing the 
central star are luminosity (L), effective temperature (Tcff), 
radius (-R,), mass (M), mass loss rate (M). 

The structure of the wind is described by the behavior 
of velocity, temperature, density, optical depth, as a function 
of the radial distance from the surface of the star. 

The equation of momentum conservation determines 
the radial variation of velocity: 



dr r^ 



(1) 



where F is the key-quantity that determines whether the 
wind is accelerated by radiation pressure. It is given by the 
ratio between the radiation pressure on the dust grains and 
the gravitational pull: 



kL 



AtvcGM 



(2) 



where k expresses the opacity, as determined by the interac- 
tion of the various species of dust grains and the electromag- 
netic radiation. In the equation describing the radial varia- 
tion of velocity we neglect pressure forces, because these are 
negligible compared to the gravity term. This assumption 
holds in the present treatment, because we do not consider 
any shock structure of the outflow. 

The description of the radial behavior of the wind is 
completed by the equation for the radial decay of the optical 
depth: 



dr 



Ri 



dr r^ 



(3) 



The boundary condition for the optical depth is that it van- 
ishes at infinity. We keep the velocity constant, and equal to 
the local sound speed, up to the region where dust begins to 
form. Whenever dust is formed, and the gas is accelerated by 
radiation pressure, the results are practically independent of 
the assumed initial value of velocity, because the asymptotic 
behavior turns out to depend only on the amount of dust 
formed. 

Finally, two additional equations describe the radial 
variation of density and temperature: 



M = 47rr pv 



Tn4 ^ rf^4 



1- 



, Ri 3 



(4) 



(5) 



The system of equations given so far is not closed, because 
the calculation of the opacity, k, requires knowledge of the 
kind of dust formed in the wind; it is therefore mandatory 
to follow the formation and growth of the dust grains of the 
various species. 
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Table 1. Dust species considered in the present analysis, their formation reaction and the corresponding key elements (see text). 



Grain Species Formation Reaction 



Key element Sticking coefficient 



Olivine 23;Mg +2{l-a:)Fe+SiO+3H20 -^ Mg2a;Fe2(i_3;)Si04 + 3H2 

Pyroxene xMg +(l-x)Fe+SiO+2H20 -^ Mg^Fe(i_^)Si03 + 2H2 

Quartz SiO + H2O ^ Si02(s) +H2 

Silicon Carbide 2Si + C2H2 -s> 2 SiC + H2 

Carbon C -5- C(s) 

Iron Fc -^ Fc(s) 



Si 


0.1 


Si 


0.1 


Si 


0.1 


Si 


1 


c 


1 


Fe 


1 



Dust growth takes place by vapour deposition on the 
surface of some pre-formed seed nuclei, assumed to be nano- 
meter sized spheres. The precise chemical nature of these 
seed grains is not important in what follows. We consider 
various types of dust, depending on the surface chemical 
composition of the star. In oxygen-rich winds, we consider 
olivine, pyroxene, quartz and iron grains, whereas for C-rich 
environments we account for the presence of solid carbon, sil- 
icon carbide and iron grains. For each condensate, we define 
a key element, whose abundance is the minimum among all 
the elements necessary to form the corresponding dust ag- 
gregate. Silicon is the key species for silicon carbide, whereas 
iron and carbon are the key elements for iron dust and solid 
carbon. As for the silicates, the key element is silicon, which 
is however replaced by oxygen when strong HBB conditions 
favor a strong depletion of the surface oxygen content. All 
the species considered, with the corresponding condensation 
reactions and key elements, are listed in Table 1. 

The growth of the size of the dust grains of a given 
species i is determined via a competition between a pro- 
duction term (Jf^), associated to the deposition of new 
molecules on already formed grains, and a destruction factor 
{Ji^'^), proportional to the vapor pressure of the key species 
on the solid state: 

^^VoAJf-Jtl (6) 

where Vo.i is the volume of the nominal molecule in the solid, 
and the growth and destruction rates are given by 



CtiUiVthA 



J"" 



OiiVth, 



'kT 



(7) 



(8) 



where rii is the number density of the key species in the wind, 
Vth,i is the thermal velocity of the same species, and Pv.i is 
the vapor pressure of the key element over the solid state. As 
we see from Eq. [SI [T] [H] the rate of growth of the dust grains 
is proportional to the sticking coefficients ai, which vary 
between zero and u nity. The coefficients ad opted for each 
species, taken from iFerrarotti &: Gaill (|2006l ). are reported 
in Table 1. 



3 THE EVOLUTION OF LOW-Z STELLAR 
MODELS 

Models with initial metallicity Z = 3 x 10~^ and mass in 
the range IMq ^ M ^ 7.5Mq are evolved through the 
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Figure 1. Variation of the temperature at the bottom of the 
convective envelope (bottom) and of the luminosity (top) of 
models of initial mass in the range 1.5Mq ^ M ^ 7.5Mq evolved 
during the AGB phase. The total mass of the star, reported in 
the abscissa, is decreasing during the evolution. 



thermal pulses phase, until most of the mass in the envelope 
is lost. We find that, for this Z, 7.5AIq is the largest mass 
not undergoing supernova explosion. 

The physical properties of the models, with the maxi- 
mum luminosity and temperature at the bottom of the en- 
velope reached during the AGB evolution, are reported in 
Table [3 

The top panel of Fig. [l] shows the variation of the lumi- 
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Table 2. Relevant properties of AGB models. The last column indicates the physical process dominating the surface chemistry. 



M/Mq 



„l (Myr) M^oic/Mq 



bcc 



(lOfiK) log(L/L0) 



Z=3 X 10"* 



1.0 


5921 


.60 


4 


3.78 


TDU 


1.5 


1065 


.68 


4 


3.98 


TDU 


2.0 


805 


.72 


11 


4.13 


TDU 


2.5 


455 


.78 


38 


4.27 


TDU 


3.0 


294 


.84 


79 


4.42 


HBB + TDU 


3.5 


207 


.87 


89 


4.54 


HBB + TDU 


4.0 


155 


.90 


98 


4.63 


HBB + TDU 


4.5 


120 


.93 


100 


4.70 


HBB + TDU 


5.0 


96.0 


.97 


111 


4.75 


HBB + TDU 


5.5 


79.0 


1.00 


117 


4.85 


HBB 


6.0 


67.0 


1.05 


121 


4.97 


HBB 


6.5 


57.0 


1.16 


124 


5.05 


HBB 


7.0 


49.5 


1.25 


127 


5.12 


HBB 


7.5 


43.5 


1.34 


141 


5.20 


HBB 



nosity of models of different mass, during the AGB phase. 
The evolution is plotted as a function of the mass (decreas- 
ing with time) , to better understand the main properties of 
the stars when most of the mass loss occurs. 

More-massive models, with larger cores, evolve to larger 
luminosities. For each mass the overall energy release first 
increases during the early AGB, then decreases in the late 
phases. While the initial increase is due to the growing core 
mass, the following drop in the luminosity is favoured by 
mass loss, which decreases the mass in th e envelope. Be - 
cause mass loss increases with luminosity (|Blockeil 119951 ). 
the models evolving at large luminosity are also those losing 
mass at the highest rates. 



3.1 The evolution of high— mass AGB stars 

The bottom panel of Fig. [T] shows the evolution of the 
temperature at the bottom of the convective zone (Tbco)- 
A clear gap separates models with mass above 2Mq from 
their smaller mass counterparts. This is the signature of 
Hot Bottom Burning, operating when the core mass ex- 
ceeds the threshold mass of ~ 0.8Mq, reached for initial 
masses M > 2.5Mq. HBB is a consequence of the steep- 
ness of the temperature gradient: under these conditions 
the bottom of the convective envelope partially overlaps 
with the CNO burning shell, thus allowing some nuclear 
activity to take place within t he outer convective zone. 
The pioneering investiga t ions b y lRenzini fc Volil (|l98lh and 
[Blocker fc Schoenberneij (|l991l ) showed that HBB ignition 
can be achieved whe never a high-efficiency treat ment of 
convection is adopted. IVentura fc D'Antonal (|2005al ) showed 
that HBB is easily obtained in models with M ^ 3Mq when 
the FST treatment of convection is used. HBB is a complex 
phenomenon, where the thermodynamic coupling between 
the internal, degenerate core and the external envelope plays 
a major role: given these arguments, HBB cannot be de- 
scribed within a synthetic modeling context, which is one of 
the main motivations of the present investigation. 

HBB can impose a strong variation in the surface chem- 
istry, because of the proton-capture nucleosynthesis experi- 
enced at the bottom of the convective zone, and the rapidity 



of convective motions that homogenize the whole envelope: 
the chemistry of stars experiencing HBB reflects the equi- 
libria of CNO nucleosynthesis. 

The variation of the surface oxygen, shown in the top 
panel of Fig. [21 is a clear indicator of the mechanism domi- 
nating the surface chemistry. Models with mass M > 5Mq 
destroy oxygen via HBB; whereas lower mass stars, with 
M < 3Mq increase it via repeated TDU episodes. Note 
that in models with mass in the range (3 — 5)Mq oxygen is 
initially destroyed by HBB, then produced by TDU episodes 
in the more advanced AGB phases (see last column of Table 
2). The oxygen content in the gas ejected by super-AGBs 
is larger than in th e matter expelled by massive AGBs: this 
efTect, discussed in lVentura fc D'Antonal l|201ll '). is a conse- 
quence of the velocity with which super-AGBs loose their 
envelope, the rate at which the mass is lost exceeding the 
rate at which oxygen is destroyed. The largest depletion of 
oxygen, in analogy with higher metallicity models (see pa- 
perl and paperH), is found for the masses at the edge be- 
tween the AGB and the super- AGB regimes, i.e. M ~ 6Mq . 

The bottom panel of Fig. [2] shows the evolution of the 
surface Si/O ratio; this will be important when evaluating 
the possibility of silicate formation in the atmospheres of 
these objects. In models with M ^ 5Mq the simultaneous 
depletion of oxy gen and production of silicon via the Mg-Al 
nucleosynthesis (jVentura. Carini fc D'Antonal 120111 ) makes 
silicon more abundant than oxygen at the surface. 



3.2 The carbon— star stage in low— mass AGBs 

HBB is not active in low-mass stars with A/ < 3Mq; the 
only mechanism able to modify their surface chemistry is 
TDU, which determines an increase in the surface carbon. 
The amount of ^^C carried to the surface after each ther- 
mal pulse depends on the extent of the inwards penetra- 
tion of the convective zone: the more internal is the re- 
gion reached by the base of the envelope, the larger will be 
the carbon increase. Unfortunately, the extent of the third 
dredge-up is sensitive to many ingredients used in calcu- 
lating the evolutionary se quences, still unknown from first 
principles (|Karakaal201in . The method used to determine 
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Figure 2. Bottom: Variation of the surface oxygen of AGB mod- 
els of various initial mass during the AGB evolution. Massive 
models, with mass M > SMq, destroy oxygen, whereas low-mass 
stars, with mass below SMq , owing to repeated third dredge-up 
episodes, increase ^^O. In models with SMq < M < &Mq oxygen 
is initially destroyed at the surface by HBB, then some produc- 
tion occurs when TDU occurs. Top: evolution of the Si/O ratio 
in the same models shown in the left panel. The dotted line indi- 
cates the stage when the number densities of the two species are 
equal. 



the location of the border between convective and radiative 
regions is crucial for the results obtained. While sticking 
to the Schwarzschild criterium underestimates the dimen- 
sions of the mixed region, the extent of the possible over- 
shoot is still an object of debate. Other mechanisms, dif- 
ferent from convective overshoot, may further extend the 
dimensions of the mixed zone: rotational mixing and inter- 
nal gravity waves may also favor mixing of layers radiatively 
stable (Stancliffc & Lattanzio 2011). 

Studies focused on the luminosity function of the car- 
bon stars in the Magellanic Clouds (IMarigo et al.l 1 1999 : 
llzzard et al.l 12004 1 calibrated the extent of TDU needed to 
reconcile the observational evidence with theoretical mod- 
eling: these investigations provided the minimum core mass 
at which TDU begins to be effective in producing carbon 
stars, and the extension of the mixed region. However, we 
are still far from a robust and reliable description, holding 
for all masses and metallicities. We decided in this context 



to avoid any detailed calibration of the possible extra mix- 
ing from convective zones present in the interiors of AGB 
stars; this analysis would be harder than for higher metal- 
licity models, from which many more observations are avail- 
able. We simply explore how the dust formation process is 
sensitive to the depth of TDU. 

The extent of TDU was changed by assuming a variable 
overshoot from the borders of the convective shell that forms 
during the thermal pulse. Accounting for this extra mixing is 
known to enhance the strength of the pulse, and to favor the 
inwar ds penetration of the c onvective envelope during the 
TDU (JHerwig fc Austin|[2004 ) . In the following, we compare 
the results obtained with no extra mixing, with those found 
when values of the overshoot parameter (see section [!])(" = 
0.001,0.002,0.005 are adopted. 

The typical evolution of a low-mass star experiencing 
TDU, in terms of the progressive increase in the surface car- 
bon, is shown in the right panel of Fig. [S] which refers to 
a 2.5Af0 model. A sudden increase in the ^^C mass frac- 
tion follows each TDU episode. The envelope of these mod- 
els becomes more and more enriched in carbon, until the 
mass left in the envelope is so small that no TDU occurs 
jVassiliadis fc Wood! [1993 '). The cases C = (black, dotted), 
C = 0.002 (red, solid) and ( = 0.005 (blue, dashed) are 
shown. Fig. U shows the average carbon abundance in the 
envelope during the AGB evolution, for models with initial 
masses IM© < M < 2.5M0. 

The average carbon in the gas ejected by the M = 
2.5M0 model is X(C) ~ 5 x 10"^, independently of (. The 
reason for this can be understood on the basis of the left 
panel of Fig. O showing the variation of the temperature 
at the bottom of the convective envelope. In the ( — 0.005 
case (blue, dashed line) the strength of the pulses is en- 
hanced such that eventually HBB conditions are reached, 
and the temperature exceeds ~ 50 x lO^K. The impact of 
( on the increase in the surface carbon is limited in mod- 
els with mass approaching the threshold for HBB ignition, 
because the increase in the strength of the thermal pulses 
eventually leads to HBB conditions, which determine a par- 
tial destruction of the carbon previously accumulated. This 
can be seen in the right panel of Fig. [3l showing that in the 
( = 0.005 model the increase in the carbon mass fraction 
due to TDU is followed by the destruction of the surface 
carbon via HBB. 



4 DUST PRODUCTION IN 

LOW-METALLICITY AGB STARS 

Dust production is strictly correlated with the surface chem- 
ical composition of the star. In our previous investigations 
we examined the metallicities Z — 0.001 (paper I) and 
Z — 0.008 (paper II): in both cases we found that in the 
winds of stars with an initial mass above the limit for the ig- 
nition of HBB, i.e. (2.5 — 3)M0, conditions are favourable for 
the formation of silicates and very little or no carbon grains 
are formed. Clearly this behaviour cannot be extended to 
extremely low metallicities, because the silicon available in 
the envelope scales with Z; we expect that below a given Z 
silicate formation is scarcely efficient, and the little amount 
of dust formed is not capable of accelerating the wind via 
radiation pressure on the dust grains. 
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Figure 3. Left: Evolution of the temperature at the bottom of the convective envelope of AGB models of 2.5M0 calculated with 
different prescriptions of the overshoot from the convective shell: no overshoot (black, dotted), C, = 0.002 (red, solid), C, = 0.005 (blue, 
dashed). The thin horizontal line indicates the limit in temperature leading to HBB conditions, with the depletion of the surface carbon. 
Right: The change in the surface carbon in the same models shown in the left panel. The effects of HBB can be seen in the decrease in 
X{^^C) during the late inter— pulse phases of the C, = 0.005 model. Open diamonds indicate the surface carbon content averaged during 
the whole AGB evolution, which corresponds to the quantity shown in Fig. (4] 



Table 3. Dust production by low-mass AGBs with Z = 3 X 10 



M/Mq 


Mc{C.=0)/Mq 


A/c(C = O.OOl)/M0 


Mc(C = O.OO2)/M0 


Mc(C = O.OO5)/M0 


Z=3 X 10-* 


1.0 
1.5 
2.0 
2.5 


1.5 X 10-* 
3.5 X 10-* 


10-* 
5.8 X 10-* 
4.5 X 10-* 


10-* 

4x 10-* 

6 X 10-* 

4.7 X 10-* 


1.5 X 10-* 
4x 10-* 

6.4 X 10-* 
5 X 10-* 



4.1 Silicate formation 

To test the possibility that sihcates are formed at the metal- 
licity investigated in this work we focus on F (Eq.[2]): the con- 
dition for the wind to be accelerated by radiation pressure 
is F > 1, which requires a minimum opacity, kmm = !klL££M.^ 

The most favorable conditions for the formation of sil- 
icates are found in the early super-AGB phases of models 
with mass M ~ 7Mq. The high luminosities reached by 
these stars (see Table [2] and the right panel of Fig.[l| favour 
large mass loss rates, which, in turn, enhance the dust for- 
mation rate, owing to the increase in the density of the wind 
(Eq. |4}. Also, as shown in Fig. [2] in these models the sur- 
face oxygen is never completely destroyed; this allows the 
formation of water molecules, present in all the reactions 
occurring when the silicates are formed (see Table 1). These 
conditions are not met during the evolution of smaller mass 
stars, whose luminosity is smaller, and that suffer a stronger 
reduction of the surface oxygen (see Fig. [2]). 

We focus on a model of 7Mq, in the phase of maxi- 



mum luminosity [L = 1.3 x lO^L©), when it has lost IMq 
(see right panel of Fig. [T|; the surface temperature is 3700 
K. Substituting M = 6Mq and L = 1.3 x 1O'^L0 in the 
expression given above we find kmin = 0.6cm'^/g. On the 
other hand, by u sing t he optical constants for silicates by 
lOssenkopf et al.l (|l992l ). in the unrealistic case that all the 
silicon is condensed into dust, we find k = 0.3cm'^/g, i.e. 
k ~ ikmin- This rcsult, compared with those from our pre- 
vious investigations, shows that Z ~ 10 is the threshold 
metallicity, below which no purely dust-driven silicate winds 
may occur. Some little production of silicates is still possible 
in these conditions, but the predictions in this case are less 
reliable, because the results depend critically on the assump- 
tions concerning the velocity with which the wind particles 
reach the condensation layer. 



4.2 Carbon and SiC dust in low— mass AGBs 

In analogy with the results for higher metallicities, in mod- 
els with mass below 3M0 dust formation depends on the 
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Figure 4. The average abundance of carbon in tiie gas ejected 
by low-mass models, whose surface chemistry is changed by the 
effects of the Third Dredge-up. The lines in the figure indicate the 
results obtained by changing the extent of the overshoot from the 
borders of the convective shell formed during the thermal pulse. 
Stellar models with M ^ 1.5M0 and no overshoot never become 
carbon stars. 



extent of TDU. The excess of carbon with respect to oxygen 
allows formation of solid carbon and silicon carbide. This 
condition is reached after a few TPs, provided that TDU is 
sufficiently deep to enrich in carbon the surface layers (see 
right panel of Fig. |3]). In the previous phases, as far as the 
C/O ratio is below unity, only silicate formation is possible. 
The percentage of silicon in the wind is however so low, that 
the rate of growth of the grains is very small, the maximum 
dimensions hardly reaching 0.05/im. In agreement with the 
previous discussion, we find that the quantity of dust formed 
is far below the threshold needed t o accelerate th e wind. This 
is consistent with the analysis by iHofnen l|2008l ). indicating 
that in the range of masses and luminosities typical of the 
low-mass regime investigated here a mimimum dimension 
of silicates grains of ^ 0.3^m is required to accelerate the 
wind. 

Table [3] reports the mass of solid carbon produced by 
AGB models of different initial mass, calculated with C, = 
0, 0.001, 0.002, 0.005. The same data are shown in Fig. 

Models with mass below 2Mq become carbon stars only 
for (^ ^ 0.002. For a given C,, the mass of carbon dust pro- 
duced increases with the mass of the star, because more 
massive models experience more mass loss and more ther- 
mal pulses, thus reaching higher carbon abundances. The 
trend in Fig. [5]reffects the behavior of the carbon content in 
the gas ejected by AGBs, shown in Fig.m The 2.5M0 model 
is found to produce a smaller amount of carbon, which can 
be explained based on the discussion in previous sections, 
concerning the HBB ignition in this model. 

Fig. [S] shows that the mass of carbon dust produced be- 
comes approximately independent of C,, for (" > 0.001. The 
main reason for this behaviour is the rapid acceleration ex- 
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Figure 5. Mass of carbon dust formed in the winds of AGB 
stars of various mass. The different symbols refer to models cal- 
culated with different extent of the extra mixing from the borders 
of the convective shell which forms during the thermal pulse. The 
shaded area indicates the degree of uncertainty associated with 
the carbon dust. 



perienced by the wind when carbon dust is produced: the 
larger is the content of carbon dust, the larger is the accel- 
eration of the wind, which prevents additional formation of 
dust. In addition, models more enriched in carbon are those 
experiencing the deepest TDU, thus evolving at smaller core 
masses: their luminosity, and consequently their mass loss 
rate, is smaller, such that the density of the gas in the wind 
is also reduced. 

Fig. |6]shows the variation of the size of C and SiC grains 
during the AGB evolution of stars of initial mass I.SMq, 
2Mq, 2.5Mq; the different lines indicate the ^ adopted in 
the computations. 

The size of the SiC grains is independent of (, because 
the key element, Si, is unaffected by TDU. Owing to the 
scarcity of silicon available, the production of SiC is modest, 
the size of SiC grains being asic < 0.02^m. The amount 
of SiC formed increases with mass, because more massive 
models experience a larger mass loss rate, that favours gas 
condensation. 

Iron is scarcely produced in all models, as a consequence 
of the small iron content of the winds of AGBs at these 
low metallicitie s. This result is consiste nt wi th previous in- 
vestiga tions by iFerrarotti fc Gaill l|2006l ) and IVentura et al.l 
(|2012ah . indicating that the iron dust produced by AGBs 
diminishes with decreasing metallicity. This picture could 
ch ange at higher metallicit ies, as suggested by the analysis 
by iMcDonald et al.l l|201lh , focused on dust production in 
stars belonging to the Globular Cluster 47Tuc. The authors 
find that iron is the main component of the dust content sur- 
rounding the two brightest sources, suggesting a favourite 
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Figure 6. The dimension of carbon and SiC dust grains formed around models with initial mass I.SMq (Left), 2Mq (Center), 2.5Mq 
(Right) during the AGB evolution. Different lines refer to the assumptions concerning the extra mixing from the borders of the convective 
shell which forms during the thermal pulses, indicated by the parameter ^ (see text for details). Dotted: no-overshoot; solid: C, = 0.002; 
dashed: C = 0.005. 



channel for iron production compared to silicates. Detailed 
models of low-mass stars at higher metallicities are required 
to be compared with these findings. 

The growth of carbon grains varies with C,, because the 
surface carbon enrichment changes according to the depth 
of TDU, as shown in Fig. |4l The grain size increases dur- 
ing the evolution, not only because more and more carbon 
is brought to the surface, but also because the luminosity, 
hence the mass loss rate, increases as the core grows in mass. 
While a minimum value of C, is required to allow TDU to 
operate efficiently, the results in terms of the size of carbon 
grains formed become insensitive to C^ as far as a sufficient 
amount of carbon, X{C) > 5 x fO"'', is present in the enve- 
lope. 

A calibration of (" is beyond the scopes of the present 
work. However, in the analysis of the carbon grain size 
formed, we may disregard the no-overshoot assumption, 
given the physical inconsistency, and the severe prob- 
lems met when comparing the results with the observa- 
tions of stars in th e Magellanic Clouds tMarigo et al. 1999; 
llzzard et al.M2004l ). Limiting our attention on the C, ^ 
cases, we conclude that the surroundings of low-mass, low- 
metallicity AGBs are expected to host carbon grains, whose 
dimensions are in the range 0.08/im < ac < 0.12/im (see 
Fig.©. These results are approximately independent of the 
initial mass of the star. 

The reason for this is the rapid acceleration of the wind 
favoured by the condensation process; the increase in the 
velocity favours a rapid decrease in the density of the gas, 
that inhibits further formation of dust. In Fig. [7] we com- 
pare the structure of the wind of two models of 2Mq , cal- 
culated with C = and C, = 0.005. The figure refers to an 
evolutionary phase when the stars have lost IMq, and the 
surface carbon abundances are Xc ~ 0.003 ((" — 0) and 
Xc = 0.014 (C = 0.005). The radial profiles of the grain 
size (red) and velocity (blue) are shown. Solid tracks refer 
to the no-overshoot models, whereas dashed lines indicate 
the C = 0.005 case. 

Dust condensation begins at a distance of ~ lOi?, from 
the stellar center, and continues until the velocity of the 



wind reaches ~ lOkm/sec. The acceleration of the wind is 
smaller in the C = case, which allows dust formation to 
occur in the expanding wind up to more external regions. 
As clearly evident in the internal panel of Fig. [T] showing 
in more details the dust-forming region, this effects partly 
compensates the difference in the carbon available in the 
winds, reducing the gap between the size of the carbon grains 
formed. 

The optical properties of solid carbon, its ability to ab- 
sorb radiation and accelerate the wind, make the description 
of the solid carbon growth only mildly dependent on the 
quantity of carbon in the stellar envelope, thus rendering 
more robust the results obtained, in terms of the dimension 
of the grains formed. 

The possibility of comparing these results with the ob- 
servations is undermined by the scarcity of data of AGB 
stars in metal poor environments. A promising field in this 
context is the increasing amount of available spectra of 
AGBs in globular clusters, from which the species and the 
quantity of dust formed can be deduced. Because masses 
evolving in Globular Clusters are surely below IM0, this 
approach offers the only opportunity to test our description 
at the lower end of the range of masses investigated, whereas 
no information can be obtained for models of higher mass. 
An example of these stud i es is th e collection of Spitzer data 
presented by ISloan et al.l (|2010h . showing a positive trend 
with Z of the amount of dust formed in the circumstellar 
shells of AGB stars. 

W hile stars in the sample presented in I Sloan et al.l 
(J2010|) are more metal rich than the models presented 
here, the two sources b elonging to NGC 4372, analyzed in 
iMcDonald et al.l l|2013l ). have a similar metallicity ([Fe/H] = 
—2.2), which offers the opportunity of a more detailed com- 
parison. From the spectra, the authors argue that the two 
stars are oxygen-rich, and dusty. The luminosity of the two 
stars is slightly below the value found in the IM0 model 
presented here (see Table 2), in agreement with the smaller 
mass (~ O.9M0) expected to evolve in the AGB of this clus- 
ter. FigO shows that some carbon-dust production is ex- 
pected in these low-mass objects, provided that some extra 
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Figure 7. Radial profiles of velocity (blue) and size of carbon 
grains (red) in the wind expanding from the surface of a 2Mq 
model calculated with no extra mixing (solid) and with f = 0.005 
(dashed). The internal panel, focused on the dust forming region, 
shows only the grains size. 



mixing from the convective zones is adopted. This is however 
no t in contrast with the fact that the two stars examined 
bv lMcDonald et al.l (120131 ) are oxygen-rich, because even in 
case that a large extra mixing is used, the model achieve the 
C-star stage only in the lates evolutionary phases, while the 
surface remains oxygen-rich for most of the duration of the 
AGB phase. 



4.3 Dust mass produced by AGB stars: the trend 
towards metal—poor chemistries 

As discussed above, the quantity of silicon present in low 
metallicity AGBs is so small to prevent the formation of 
silicates, if not in extremely small quantities. Even in the 
super-AGB regime, where the strong mass loss might favor 
the condensation process, the density of silicon in the wind is 
too small to allow dust production. Our AGB mod els expe- 
rienc e stronger HBB than other literature results (JKarakag 
uOld ). and are thus those most favorable to the formation of 
silicates. This makes these results more general, allowing us 
to conclude that Z ~ 10~^ is the threshold limit in metal- 
licity, below which no meaningful silicate formation occurs 
in the winds of AGBs. 

The shortage of silicon in the wind also affects the for- 
mation of SiC, whose production rate scales with the silicon 
density. 

The comparison of the results obtained in the present 
investigations with those presented in paper I and paper II 
allows us to infer how the mass and composition of dust 
produced by AGBs change with metallicity. Fig. [8] shows 
the mass of dust produced by stars with masses in the range 



Figure 8. The mass of dust produced by AGB models of dif- 
ferent initial mass, at various metallicities. Full points indicate 



results for Z 



IQ- 



published in paper I, open squares refer 



to the Z = 8 X 10~ metallicity presented in paper II, whereas 
full triangles refer to the present investigation. The models in the 
low-mass regime were calculated by assuming an extra mixing 
f = 0.002 from the convective shell which forms during each TP. 
The thin, dotted line separates models producing silicates from 
those producing carbon dust. 



IM© ^ M < SM© and metallicities Z = 10^^ (paper I), 
Z = 8 X 10"^ (paper II), and 3 x 10"** (present study). 

Stars with mass above ^ 3Mq produce silicates, in 
quantities scaling with Z. Low-mass stars produce mainly 
carbon dust, although some silicate formation occurs, in 
quantities increasing with metallicity, owing to the larger 
silicon available. As far as the C-star stage is reached, the 
carbon-dust formed exceeds by far the previous production 
of silicates, because of the large excess of carbon with re- 
spect to oxygen favoured by repeated TDU episodes. Most 
of the mass of dust shown in Fig. |S1 produced in the low- 
mass regime, is under the form of carbon-dust. 

The effect of metallicity on the carbon formed is not 
trivial, because the carbon dredged to the surface during 
the TDU is produced in the He-burning shell, and is thus 
independent of Z. In \ow-Z m odels TDU is more efficient 
(jBoothrovd fc Sackmanrj|l988l ). thus their surface carbon 
is larger. This effect is compensated by the cooler surface 
temperatures at which higher-Z models evolve: the dust- 
formation zone is closer to the stellar surface, thus the con- 
densation process is more efficient. Moving to smaller metal- 
licities the range of masses involved in the solid carbon for- 
mation becomes progressively smaller, because the larger 
core masses favour HBB conditions for smaller and smaller 
initial masses. For masses M < 2Mq, lower Z models pro- 
duce more dust, because they reach earlier in their evolution 
the carbon star stage. 

The results shown in Fig. [8] indicate that while the pro- 
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duction of silicates turns out to be extremely dependent on 
metallicity, the carbon-dust produced show only (if any) a 
modest dependence on Z, but the lower and upper limits of 
the masses involved. The se findings are in agreement with 
the results presented by ISloan et alj l|2012l) , showing that 
that the amount of amorphous carbon produced is fairly in- 
dependent of metallicity, for carbon stars observed not only 
in the Galaxy and in the Magellanic Clouds, but also in 
other galaxies of the Local Group. 

We plan to further extend these results to smaller metal- 
licities, but the general pattern can be inferred from the re- 
sults found here and from the properties of AGB stars at 
extremely small Z. In the high-mass domain, silicate for- 
mation becomes harder with decreasing Z, because of the 
smaller amount of silicon available and the stronger HBB, 
which triggers the destruction of surface oxygen. For masses 
M > 3Mq HBB is expected to occur down to Z = 10~* 
(|Lau et al.ll2009l V which prevents the formation of carbon- 
type dust. 

At Z below 10"*^, low-mass AGBs evo lve differently 
JLau et al.ll2009l : ICampbell fc Lattanzidl2008l '). The ignition 
of He-burning, both in the core or in a shell, is accompanied 
by the development of an extended convective zone, which, 
owing to the small height of the entropy barrier, penetrates 
the H-rich region. The following inwards penetration of the 
external envelope allows mixing of processed material in the 
interior with the surface, which is thus enriched in CNO ele- 
ments. The composition of the winds of these stars depends 
on which mechanism is more efficient in changing the surface 
chemistry. Stars with mass above ~ IMq experience both 
HBB and TDU; although dual shell flash is also likely to 
occur, the change in the surface chemistry is mostly deter- 
mined by the combined effects of TDU and HBB: the former 
favors an enrichment in the surface carbon after each TP, 
the latter destroys part of this carbon and produces nitro- 
gen. As a result, the average carbon mass fraction barely 
reaches X{C) = 5 x 10~*, thus preventing the formation of 
meaningful quantities of carbon grains. Therefore, we can 
expect that when Z ^ 10~^, dust production is limited to 
low-mass stars, whose evolution is not influenced by HBB 
and whose surface carbon may survive to proton fusion. 



4.4 The uncertainties in the modelling 

Before addressing the implications of the results presented, 
we discuss the robustness of these findings. The treatment 
of convection and mass loss are known to have a significant 
impact on AGB modelling, both on th e chemical and phys- 
ical s ides I Ventura fc D'Antonall2005al lbl: [Ventura fc Marigol 
I2OI0I ). However, at odd s with previous explo rations based 
on higher metallicities (jVentura et al.l l2012al la) , dust pro- 
duction in the high-mass domain is not sensitive to these 
uncertainties at these low Z's: the amount of silicon in the 
wind is so scarce, that no meaningful silicate-type dust is 
expected, far below the quantities needed to accelerate the 
wind. 

The choice of the convection modelling influences the 
upper limit of the masses producing carbon-type dust: if 
a less efficient description of convection was adopted, the 
threshold limit (in mass) to experience HBB would shift 
~ O.5M0 upwards, which would extend the range of masses 
involved in carbon production up to ~ SM©. 



The description of mass loss could potentially influence 
the physical evolution of low-mass stars: the mass loss rate 
determines the number of thermal pulses experienced, thus 
the extent of the surface carbon enrichment. Yet, we do not 
expect a strong sensitivity of the amount of carbon-dust 
formed to the mass loss description for the following rea- 
sons: a) a higher carbon enrichment favours a rapid cooling 
and expansion of the outer layers of the stars, with the con- 
seque nt increase in the mass loss rate experienced ([Marigol 
[2003), which prevents many more TDU episodes to occur; b) 
as discussed previously, the size of the carbon grains formed 
are not strongly sensitive to the amount of carbon present in 
the wind, owing to the high opacity of carbon grains, that 
favour a strong acceleration of the winds as soon as they 
begin to form. 

Concerning the wind modeling, the robustness of the 
present results must be confronted with a more complete 
description, which accounts for the presence of shock waves. 
Dust formation would be favoured by the shock, because 
compression of the gas lifted upwards would increase the 
densities of the key species, and dense parcels of gas would be 
pushed to cool regions distant from the central star, where 
dust is expected to be formed in great quantities. 

In oxygen-rich environments, a further effect due to 
the presence of shock waves might be the deviation from 
chemical equilibrium in the gas phase, with the consequence 
that carbon is not entirely locked into CO molecules, part 
of it remaining available for carbon-dust formation. This 
scenario could lead to the simultaneous formation of carbon 
and silicates grains in the circumstellar shells of oxygen- 
rich AGBs, and would provide an alternative mechanism 
for a dust driven mass loss for M stars tHofner fc AndersenI 
l2007f) . 



5 CONCLUSIONS AND IMPLICATIONS FOR 
COSMIC DUST ENRICHMENT 

We calculate the dust formed around AGB and super-AGB 
stars of metallicity Z = "i x 10~*, the lowest metallicity for 
which dust production by low and intermediate mass stars 
has been considered so far. This study complements previous 
investigations that were limited to metallicities Z — 10~^ 
and Z — 9, X 1Q~'^ . Stellar evolution is followed by means 
of a full integration of the equations of stellar structure, al- 
lowing to describe self-consistently the Hot Bottom Burning 
phenomenon. 

We find that for stars with initial metallicity Z = 3 x 
10"*: 

• In the high-mass domain, with M ^ ZMq, strong HBB 
prevents the formation of carbon stars. The scarcity of sil- 
icon in the wind and the strong destruction of the surface 
oxygen in the most massive models lead to a negligible for- 
mation of silicates, not sufficient to accelerate the wind via 
radiation pressure on dust grains. 

• Dust production is limited to the formation of solid 
carbon grains around low-mass stars, with initial mass 
M ^ 2.5M0. The results depend on the way TDU is mod- 
elled, but as far as a minimum abundance of carbon of the 
order of Xc ^ 0.005 is reached in the surface layers, the re- 
sults are approximately independent of the extent of TDU. 
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This is a consequence of the balance between larger car- 
bon abundances reached by models experiencing the deepest 
dredge-up and their lower core masses, that lead to smaller 
mass loss rate. 

• Under these conditions, the carbon dust produced is 
2 X 1O"*M0 < M < 6 X 1O~''M0 with grains radii in the 
range 0.08/im < ac < 0.12p m. 

• Compared to higher metallicity models, low-mass stars 
with M ~ 1 — I.5M0 produce more dust because they 
reach more easily the carbon-star stage. Conversely, for 
M > I.5M0 higher Z models produce dust more efficiently, 
owing to their cooler surface layers. 

• These results can be extrapolated to even more metal- 
poor AGBs. No dust is expected in the high-mass domain, 
dominated by HBB, because of the scarcity of silicon avail- 
able. At low Z's HBB is present even in low-mass stars, 
below 2M0. Although it is not sufficiently strong to burn 
oxygen, still it destroys part of the carbon accumulated by 
TDU. The scarcity of carbon in the envelope prevents the 
formation of carbon grains, if not in small quantities. 

These results have interesting implications for the con- 
tribution of AGB stars to dust enrichment. The grid of stel- 
lar models that we have presented in paper I, II and in the 
present study, suggest that for Z = 8 x 10"'^ = OAZq, 
the metallicity of the Large Magellanic Cloud (LMC), AGB 
stars can contribute to dust enrichment of the ISM on rela- 
tively short time-scales, « 40 Myr, comparable to the evolu- 
tionary time of a 8M0 star. According to our results, the en- 
richment is initially limited to silicate dust, and carbon dust 
is released on longer timescales, t > 300 Myr, comparable to 
the evolutionary time of a 3Mq star. For initial metallicities 
in the range 3 x lO""* < Z ^ 10"^ or O.OISZq < Z ^ O.OSZq , 
the production of silicates by the most massive AGB stars 
is strongly suppressed and low-mass stars can mostly con- 
tribute to carbon dust enrichment when t > 300 Myr. 

It is important to stress that these results and the de- 
pendence of the predicted AGB and super-AGB dust yields 
on the initial stellar metallicity reflect the complex evolu- 
tion of these stars and can not be appreciated when syn- 
thetic stellar models are adopted. Indeed, our results partly 
contradict previous claims that AGB stars can always con- 
tribute to carbon dust enrichment, independently of their 
initial mass and metallicity, at least in the metallicity range, 
5 X 10'^Zq ^ Z ^ IZq (Ferrarotti & Gail 2006; Zhukovska 
et al. 2008). 

More importantly, our study suggests that the dust 
yields computed for stars with initial metallicity of Z ~ 
10~^Zq can not be extrapolated to lower metallicity as 
the contribution of AGB stars to dust enrichment can be 
safely neglected when the metallicity of the stars is Z < 
10~'*(5 X 10~^Zq). Hence, our results imply that at these 
low metallicities, supernovae are left as the only viable stel- 
lar dust sources (Todini & Ferrara 2001; Nozawa et al. 2003, 
2007; Bianchi & Schneider 2007; Cherchneff & Dwek 2010). 
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